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The steepening (break) of the power-law fall-off observed in the optical emission 
of some GRB afterglows at epoch ~ 1 day is often attributed to a collimated outflow 
(jet), undergoing lateral spreading. Wider opening GRB ejecta with a non-uniform 
energy angular distribution (structured outflows) or the cessation of energy injection 
in the afterglow can also yield light-curve breaks. 
j>^, ■ We determine the optical and X-ray light-curve decay indices and spectral energy dis- 

\ tribution slopes for 10 GRB afterglows with optical light-curve breaks (980519, 990123, 

IT) ! 990510, 991216, 000301, 000926, 010222, 011211, 020813, 030226), and use these prop- 

\»Q ■ erties to test the above models for light-curve steepening. It is found that the optical 

<^> breaks of six of these afterglows can be accommodated by either energy injection or 

l/"~) ■ by structured outflows. In the refreshed shock model, a wind-like stratification of the 

circumburst medium (as expected for massive stars as GRB progenitors) is slightly 
^ c"| favoured. A spreading jet interacting with a homogeneous circumburst medium is re- 

P ,■ quired by the afterglows 990510, 000301, 011211, and 030226. The optical pre- and 

post-break decays of these four afterglows are incompatible with a wind-like medium. 
The current sample of 10 afterglows with breaks suggests that the distribution of the 
break magnitude Aa (defined as the increase of the afterglow decay exponent) is bi- 
| modal, with a gap at Aa ~ 1. If true, this bimodality favours the structured outflow 

model, while the gap location indicates a homogeneous circumburst environment. 



Key words: gamma-rays: bursts - ISM: jets and outflows - radiation mechanisms: 
non-thermal - shock waves 



1 INTRODUCTION 

The optical light-curves of most afterglows monitored over 
more than one decade in time exhibit a steepening at about 
1 day after the GRB. Such a break has been predicted by 
the sideways spreading jet model of Rhoads (1999) and was 
observed shortly afterward for the first time in the afterglow 
990123 (Kulkarni et al. 1999). Since then, many other after- 
glows displayed an optical light-curve break, to which the 
jet model has been extensively applied. 

Rees & Meszaros (1998) have proposed that the blast 
wave may be "refreshed" when there is a substantial energy 
input in the forward shock. Such an energy injection could 
occur either if there is a wide distribution of Lorentz factors 
in the GRB explosion, the slower ejecta catching up with 
those which were initially faster during their deceleration 
by the circumburst medium, or if the central energy is long 
lived, releasing a relativistic outflow for hours or days after 
the explosion. The increase of the ejecta energy mitigates 



the deceleration of the forward shock, which results in a 
slower decay of the afterglow flux, as was proposed by Fox 
et al. (2003) to explain the nearly flat optical light-curve of 
the afterglow 021004 at 0.003-0.1 days after the burst. In 
this model, a steepening of the afterglow decay would occur 
when the rate of energy injection decreases. 

Another possibility, proposed by Meszaros, Rees & Wi- 
jers (1998), is that the angular distribution of the ejecta 
kinetic energy is not uniform. Then, as the fireball decel- 
erates, the relativistic beaming of the radiation allows the 
observer to receive emission from an ever wider region of the 
fireball, whose energy varies with observer time differently 
than in the case of an isotropic energy distribution. In this 
model, a light-curve break would result if the second order 
angular derivative of the energy per solid angle is negative, 
as in the scenario proposed by Rossi, Lazzati & Rees (2002) 
and Zhang & Meszaros (2002). 

The purpose of this paper is to test the above possi- 
ble models for afterglow beaks: refreshed shocks, structured 
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outflows, and spreading jets. To this end, we compare the 
light-curve breaks expected in each model with those mea- 
sured in the optical emission of ten GRB afterglows, and 
add consistency checks with the X-ray light-curves and the 
optical-to-X-ray spectral energy distribution (SED) slope. 



2 OPTICAL CONTINUUM SLOPE AND 
LIGHT-CURVE INDICES 

The simplest and most widely used test of the blast wave 
model for GRB afterglows consists of comparing the light- 
curve decay index with the SED slope. It is a robust pre- 
diction of the standard blast wave model (Meszaros & 
Rees 1997), confirmed by observations (e.g. Wijers, Rees 
& Meszaros 1997), that the afterglow light-curve at a fre- 
quency v above the peak of the afterglow spectrum decays 
as a power-law 



F v oc t 



(1) 



where a depends only on the location of the cooling fre- 
quency (y c ) relative to v, the slope p of the power-law distri- 
bution dN/de with the energy e of the electrons accelerated 
by the blast-wave 

dN/de oc e~ p , (2) 
and the radial stratification of the circumburst medium: 



n(r) oc r 



(3) 



with s — for a homogeneous environment and s = 2 for 
the free wind of a massive GRB progenitor. For a fireball 
(or a jet at times when it is sufficiently relativistic that rela- 
tivistic beaming prevents its finite angular opening to affect 
the observed afterglow emission), the light-curve index of 
equation (1) is given by 



where the order in the rhs is the same as in equation (4). 
Equation (7) is often used to test the blast wave model and 
to identify the type of external medium (if v < v c ). However, 
this exercise may not always provide a conclusive result be- 
cause a small amount of dust extinction in the host galaxy 
can alter significantly the slope of the optical SED. For an 
extinction of At z+1 - )1/ magnitudes at an observer frequency 
v, the local slope of the reddened SED is 



d 



d In v 



]n[F„ ■ 10" 



= /3 o + 0.9A 



din A v 
d\nv 



(8) 



where f3 is the intrinsic (unreddened) SED slope. For a 
simple A v oc v reddening curve, 



Piy) = f3 + 0.9 (z + 1) A v 



(9) 



where A v is the extinction in the host frame at the observing 
frequency v. Hence, a host extinction of Ay — 0.1 mag is 
enough to confuse the test based on equation (7) for an 
afterglow at redshift z > 1. 

There are a few mechanisms (described below) which 
can alter the light-curve decay index a of equation (4) and, 
hence, modify the test equation (7). 



2.1 Energy Injection 

As proposed by Rees & Meszaros (1998), the blast wave 
may be "refreshed" by some delayed (or lagged) ejecta. For 
convenience, we quantify the energy injection by expressing 
the fireball energy £ as a function of the observer time t. 
That afterglow light-curves are power-laws in observer time 
indicates that E(t) is not far from a power-law either. For 
this reason, we consider that 




v < v a & s — 
v c <v & s = 0, 2 
v < v c & s — 2 



(4) 



(Sari, Piran & Piran 1998). The second line in the rhs of 
equation (4) assumes that electron cooling is dominated by 
synchrotron losses. If radiative losses occur mainly through 
inverse Compton scatterings, then 



a{v > v c ) = -p- 



f(4-p)- 1 
1 

p/(8 - 2p) 
1.5 



s = 

s = 

s = 2 

s = 2 & 3 < p 



2<p<3 
3<p 
2 <p< 3 



(5) 



(Panaitescu & Kumar 2001). We ignore the possibility of in- 
verse Compton dominated electron cooling, keeping in mind 
that, because of this deficiency of our treatment, we may 
miss some viable models with v c < v. 

For the electron distribution given in equation (2), the 
intrinsic afterglow SED at optical and X-ray frequencies is 
a power-law: 



2 1 V , v c <v 



(6) 



Eliminating the electron index p between equations (4) and 
(6), leads to 



3/3 -2a = 0,1,-1 



(7) 



E(t) oc t e 



(10) 



Assuming negligible radiative losses after about 0.1 day, 
when afterglow observations are usually made, the fireball 
Lorentz factor is given by F 2 M(r) oc E, where M(r) oc r 3 ~ s 
is the mass of the swept-up circumburst medium. Together 
with r oc F 2 t, equation (10) is equivalent with the following 
distribution of the ejecta energy with Lorentz factor: 



dE 



7 : 



(7e + 3)/(3-e) , s = 
(3e + l)/(l-e), s = 2 



(11) 



for e < 3 (e < 1) for a homogeneous (wind-like) medium. 
The effect of energy injection on the afterglow flux F v de- 
cay can be easily assessed from the expression of F v (t) as 
function of the fireball energy E for the no-injection case. 
From equations (B7), (B8), and (C5) of Panaitescu & Ku- 
mar (2000), the effect of energy injection is a reduction of 
the light-curve decay index given in equation (4) by 




h>o < Vc & S = 
V c < Vo & s = 0, 2 
v < v c & s = 2 



(12) 



When the injection of energy subsides, the power-law after- 
glow decay steepens from t~( Q ~ Ac "«) to t~ a , assuming that 
the ejecta are spherical or sufficiently relativistic that the 
effects associated with collimation are not yet manifested. 
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2.2 Structured Outflows 

Another possibility is that the ejecta energy per solid angle, 
dE/dQ, is not constant, but varies with the angle 9 measured 
from the outflow symmetry axis (Meszaros, Rees & Wijers 
f998). The acceleration of the GRB outflow and penetration 
through what is still left of the progenitor stellar envelope 
should lead to a dE/dQ decreasing with increasing polar 
angle 6, as illustrated in fig. ff of MacFadyen, Woosley & 
Heger (2001), given that the mass per solid angle of the 
helium core increases with 9 (hence a higher binding energy 
and a larger dissipation of the jet energy, leading to a faster 
jet spread, at larger angles) . That afterglow light-curves are 
power- laws in time suggests that dE/dQ can be sufficiently 
well approximated as a power-law in 9: 



dE/dQ oc 



const 

(o/Oc)- q 



< ■ 
> ■ 



(13) 



with a uniform core of opening 9 C to avoid the on-axis diver- 
gence. To make use of analytical results for the light-curve 
decay index (Panaitescu & Kumar 2003) , we assume in this 
work that the direction toward the observer lies within the 
core opening. In this case, when the blast-wave has decel- 
erated enough that the emission from the outflow envelope 
is beamed relativistically toward the observer (i.e. when the 
fireball Lorentz factor F drops below 9~ x ), the afterglow 
emission decay steepens from that given in equation (4) to 
F v oc t- ( - a+Aa "'\ where 



3g(8 - g) 
45(8 - q) 

g(4 - q) 

29(4 - q) 



(14) 



because the envelope has an energy density smaller than 
that of the core. Equation (14) holds provided that q < q, 
with 



8 



(S. 
(S. 



0, v < v c ) or (s 
0, v c < v) or (s 



2, v c < v) 
2, v < v c ) 



(15) 



If the observer lies outside the core, the pre-break decay in- 
dex cannot be obtained analytically. Obviously, as the more 
energetic core becomes visible to the observer, the afterglow 
decay should be shallower than that given in equation (4) 
and the break magnitude larger than that of equation (14) 
(see fig. 2 and 3 in Panaitescu & Kumar 2003). 

For q > q the ejecta energy density dE/dQ. falls-off away 
from the axis sufficiently fast that the afterglow emission is 
dominated by the core. In this case, when the core boundary 
becomes visible the afterglow decay index increases by 



Aa j/o = 



3/4, 
1/2, 



(16) 



independent of frequency (Panaitescu, Meszaros & Rees 
1998). The subscript "j/o" stands for "jet in an outflow" 
and hereafter will identify the model of a jet whose lateral 
spreading is prevented by the envelope but which gives es- 
sentially all the afterglow emission. 

2.3 Jets 

An even larger break Aa than given in equation (16) is ob- 
tained if the envelope has too little ejecta to impede the lat- 



eral spreading of the core, as the sideways expansion of the 
jet leads to a faster decrease of its Lorentz factor, starting at 
about the same time when the jet edge becomes visible. This 
model represents the extreme case of a structured outflow 
with dE/dQ a "top-hat" function (i.e. g — > oo in eq.[13]). 
When the jet edge becomes visible, the light-curve power- 
law decay steepens from that in equation (4) to F v oc t~ p 
(Rhoads 1999), i.e. the light-curve index increases by 




v < v c & s = 
v c < v & s = 0, 2 
v < v c & s = 2 



(17) 



3 MODEL SELECTION CRITERIA 

To summarize, in the framework of the above models, af- 
terglows light-curves should exhibit a steepening from F v oc 
t- ai to F v oc t- a2 , with: 



Energy Injection (EI) 
Structured Outflow (SO) 
Jet in Outflow (J/O) 
Jet (J) 



J cti — a 
I «2 = a 



— a — Aa e 



cti = a 

c*2 = a + Aa s 



ai = a 

Q2 = a + Aaj/ 



cti = a 

<*2 = a + Aatjet 



(18) 
(19) 
(20) 
(21) 



where a is that given in equation (4). 

Equation (6) for the intrinsic SED slope f3 , equation (9) 
for the observed slope f3, and equations (18)-(21), provide a 
simple criterion to determine from the observables a and (3 
which case is at work, for a given model and afterglow: 



i) 3/3-2qi, 2 <-1 

ii) 3/3-2ai, 2 € (-1,0) 
Hi) 3/3 - 2ai, 2 € (0, 1) 
iv) 3/3 - 2ai, 2 > 1 



model does not work 
v < v c & s = 2 
v a < v c h s — 0, 2 
no restriction 



,(22) 



where v ~ 3 ■ 10 14 Hz. In the above criterion, the pre-break 
index ct\ applies to the SO, J/O and J models, and the 
post-break index e*2 is to be used for the EI model. The first 
line of equation (22) states that even the largest possible 
intrinsic SED slope, /3 = /3, is too small to be consistent 
with the steepness of the post-break light-curve decay. The 
following lines specify which cases are allowed; in each case 
the difference /3 — j3 giving the host extinction Av ■ 

Once the working cases for the EI model are identified, 
the break magnitude Aa can be used to determine the ex- 
ponent e of the injection law. From equations (4) and (12): 



e = 3 Aa 



(a 2 +3)- 1 
(aa + 2)- 1 
(aa + l)- 1 



v < v c h s = 
v c < v & s = 0, 2 
v < v c & s = 2 



(23) 



The only restriction here is e < 3 for a homogeneous medium 
and e < 1 for a wind-like medium, otherwise the dynamics of 
the refreshed shocks would be inconsistent with the assumed 
injection law E oc t e . 

Taking the break magnitude given in equation (17) as 
an upper limit for the Aa allowed by the SO, J/O and J 
models, and keeping in mind that larger breaks could result 
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icr 4 1(T 3 icr 2 io~ 1 10° io 1 10 2 io 3 10 4 

time (day) 



Figure 1. Optical (upper panel) and X-ray (lower panel) light-curves for 10 GRB afterglows with breaks and power-law fits to the pre- 
and post-break emission. In most cases, the decay indices given in Table 1 are averages of the fitting indices over a few bands, but only 
one band (R) is shown here. The X-ray light-curves were calculated from the available wide band (0.2-10 keV or 2-10 keV, usually) 
fluxes and SED slopes. For clarity, light-curves have been shifted horizontally (to the left for "«— ", to the right for "— »") by the indicated 
factors. For each afterglow, the shifting factor is the same for both the optical and X-ray light-curves. The data shown here and the 
optical measurements of Fig. 3-11 are taken from : 980519 - Halpcrn ct al. (1999), Vrba ct al. (2000), Jaunscn ct al. (2001), Nicastro 
et al. (1999); 990123 - Galama et al. (1999), Kulkarni et al. (1999), Piro (2000); 990510 - Harrison et al. (1999), Stanek et al. (1999), 
Kuulkers et al. (2000); 991216 - Garnavich ct al. (2000), Halpern et al. (2000), Corbet & Smith (1999), Piro ct al. (1999), Takeshima 
ct al. (1999); 000301 - Jensen et al. (2001), Rhoads & Fruchter (2001); 000926 - Fynbo ct al. (2001), Harrison ct al. (2001), Price et al. 
(2001), Piro et al. (2001); 010222 - Masetti et al. (2001), Stanek et al. (2001), in't Zand ct al. (2001); 011211 - Holland ct al. (2002), 
Jakobsson et al. (2003), Borozdin & Trudolyubov (2003); 020813 - Covino ct al. (2003), Gorosabel et al. (2004), Butler ct al. (2003); 
030226 - Klose et al. (2004), Pandey et al. (2004). 



from structured outflows if the observer lies outside the core 
opening, one reaches the following criterion for identifying 
the viable cases for structured outflows and jets: 



i) 3 Aa — a?i > 3 
it) 3Aa - qi e (2,3) 
Hi) 3Aa - qi e (1,2) 
iv) 3 Aa — ai < 1 



observer outside core 
v < v c & s = 
s = or v c < v & s — 2 
no restriction 



(24) 



For the SO model, if the observer is located within the open- 
ing of the uniform core, the measured Aa determines the 
exponent q of the angular distribution of the ejecta kinetic 
energy (eq.[13]). From equations (4), (14) and (15), one ob- 
tains 



q — Aa 



8/(q 2 +3) 
8/(q 2 +2) 
12/(3a 2 - 2ai + 1) 
12/(3a 2 -2ai +2) 



V < V c & S — 

v c < v a & s = 
v < v c & s = 2 
v c < v & s = 2 



.(25) 



If the direction toward the observer is outside the outflow 
core, then the criterion given in equation (24) does not hold. 
From the maximum post-break decay index allowed by the 
SO model 



Ct2 = 



1 J 3p , (s = 0, v < v c ) or (s = 2, v c < v a ) 



4 I 3p + 1 , (s — 0, v c < v ) or (s = 2, v a < v c ) 



(26) 



and from (3 > 
replaced by 



the criterion given in equation (22) is 
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i) 3/3 


- 2a 2 < - 


2 




SO model does not work 


ii) 3/3 


- 2a 2 G (- 


-2,-1) 




v < v c & s = 2 


Hi) 3/3 


- 2a 2 G (- 


-1,-0.5) 




< ^ c & s = 0, 2 


iv) 3/3 


- 2a 2 G (- 


-1/2,0) 




s = or v < !/ c & s = 2 


iv) 3/3 


- 2a 2 > 






no restriction 



4 OBSERVATIONS 

Before proceeding with testing the models presented in §2, 
we determine the relevant observables for a sample of 10 
GRB afterglows whose optical light-curves exhibit a break. 
The pre- and post-break decay indices, a\ and a<i, of the 
optical light-curves were obtained through power-law fits to 
the afterglow light-curve at times before and after the break, 
excluding measurements close to the break time. Therefore, 
the fitted data stretch a time range somewhat shorter than 
given in columns 2 and 3 of Table 1. Such fits are shown 
in only one band, for each afterglow, in the upper panel of 
Fig. 1. With the exception of 000301, all pre-break indices 
Qi are averages of values determined for 3 or 4 optical bands 
(7, 7?, V, and B). However, for half of the afterglows, suf- 
ficient data to determine the post-break index exist in only 
one band (usually the R band). The lower panel of Fig. 1 
shows the power-law fits to the X-ray emission, used to de- 
termine the pre- or post-break decay indices S\ and 62 ■ With 
the exception of 010222, the existing data allows the deter- 
mination of only one of the indices; for 010222, the existence 
of a break relies on single, late-time measurement. 

Fig. 2 shows the distribution of the break magnitude 
Aa. With the current sample of ten afterglows with breaks, 
there is only a vague suggestion for a bimodal distribu- 
tion, with a gap at about Aa = 1. Such a bimodality, if 
real, would not be a natural feature of the Energy Injection 
model, where any value of Aa is (a priori) equally likely, 
but could arise from Structured Outflows, where the smaller 
break magnitudes are due to the ejecta angular structure, 
while the larger values are for spreading jets. The gap at 
Aa = 1 would then point to a homogeneous medium and 
cooling frequency above the optical. 

Figs. 3-12 show the slope /3 of the optical SED for each 
afterglow, at a few epochs. With the exception of 011211, 
these epochs bracket the light-curve break time, extending 
over approximately one decade in time. We do not find any 
compelling evidence for a spectral evolution in any after- 
glow, though, given the uncertainty of the instantaneous 
SED slopes, the crossing of a slowly evolving, shallow, and 
smooth spectral break (such as the cooling break) cannot be 
ruled out. For this reason, we do not consider in this work 
light-curve breaks arising from the passage of a spectral 
break through the observing band. Such a scenario would 
require a break in the electron energy distribution which is 
stronger than that due to cooling or which evolves faster. 
The evolution of the cooling break (y c tx i -1 / 2 for s = and 
v c oc t 1 ' 2 for s = 2, assuming synchrotron-dominated radia- 
tive losses) and the steepening A/3 = 1/2 of the SED across 
it imply that the passage of v c through the optical domain 
would yield only a modest break magnitude: Aa = 1/4. 

The optical slopes /3 reported in Table 1 are averages 
of the instantaneous values shown in Fig. 3-12. Besides re- 
ducing the uncertainty, averaging over a few epochs is also 
a safeguard against short-term colour fluctuations. Optical 



fluxes were corrected for Galactic extinction but not for a 
possible host reddening. While the curvature of the opti- 
cal SED due to the differential dust extinction produced by 
Av ~ 0.1 mag is usually masked by the errors of measuring 
the afterglow I, R, V, B band fluxes, near-infrared (K, H, J 
bands) measurements widen enough the frequency baseline 
that they allow the measurement of that curvature, the de- 
termination of Av and of /3 (e.g. Fynbo et al. 2001, Jensen 
et al. 2001, Jakobsson et al. 2003). Aside for the uncertainty 
in the fitted /3 introduced by assuming a reddening curve, a 
second problem with inferring /3 from the observed curva- 
ture of the afterglow optical SED is that the intrinsic spec- 
trum may be curved if the cooling frequency v c is within 
the NIR-optical domain, a coincidence which may not be 
that rare. In such a case, the host extinction will be overes- 
timated. 

An other way of inferring Av and /3 is to compare the 
X-ray flux with the extrapolation to the X-ray domain of 
the optical SED. If the measured X-ray flux lies on that 
extrapolation, then it is tempting to conclude that Av ^ 
0, while an X-ray excess may be seen as the signature of 
optical host extinction. In the latter case, from the optical- 
to- X-ray flux ratio and the optical and X-ray SED slopes 
(/3 and f3 x ), it is possible to determine the v Cl f3 , and Av 
whenever f3 x is known sufficiently well (e.g. Galama & Wijers 
2001). Otherwise, if v c is assumed not to be in between the 
optical and X-ray domains, matching the extrapolation of 
the dereddened optical continuum with the X-ray flux could 
lead to an underestimation of Av and overestimation of /3 . 
In the worst case - v c J; v - the X-ray flux would lie on 
the extrapolation of the optical SED for Av = 0.5/(2 + 1). 

Perhaps the largest uncertainty in determining v c , /3 , 
and Av as described above arises from the possible contribu- 
tion of inverse Compton scatterings to the X-ray afterglow 
emission, which alters the optical-to- X-ray flux ratio and 
the X-ray SED slope. If inverse Compton scatterings are 
ignored, it could lead to an overestimated Av and under- 
estimated /3 . For this reason, we do not determine /3 and 
Av directly from observations, but infer them from the ob- 
served SED slopes and light-curve decay indices, using the 
relationship between them that is expected for each model. 
For further use, we also measure and list in Table 1 the 
observed optical-to-X-ray slope, defined by: 

A- = -^, (28) 

where the subscripts "x" and "o" denote X-ray and optical. 
As we shall see, the inferred host extinction Av is usually 
less than 0.2 mag, i.e. the intrinsic optical-to-X-ray slope is 
only slightly larger than the observed one. 

The la dispersion of the indices ai, 0:2, /3, 5\, 82, and 
Pox given in Table 1 is for the joint variation of both the 
normalization constant and the index, i.e. it is the width of 
the projection of the A\ 2 = 2.30 contour on the a-axis. 



5 MODEL TESTS AND RESULTS 

We use the criteria given in equations (22) and (24) to se- 
lect from the four models for light-curve breaks - Energy 
Injection, Structured Outflow, (non-spreading) Jet in Out- 
flow, and (spreading) Jet - those which can accommodate 



© 2005 RAS, MNRAS 000, 000-000 



6 A. Panaitescu 



Table 1. Light-curve decay indices and slopes of the optical and X-ray emission of 10 GRB afterglows with optical light-curve breaks 



GRB 


ti/t b 


ti/tb 




CK2 


P 


Si 


&2 


Px 


Pox 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


980519 


0.3 


3 


1.87 ±0.11 


2.43 ± 0.35 


1.07 ±0.12 


2.23 ±0.91 




1.8 ±0.6^ 


1.01 ±0.03 


990123 


0.1 


£ 10 


1.17 ±0.05 


1.65 ± 0.20 


0.63 ±0.05 


1.44 ±0.13 






0.67 ±0.02 


990510 


0.1 


20 


0.92 ±0.04 


2.21 ± 0.09 


0.60 ±0.10 


1.40 ±0.19 




1.03 ± 0.08^ 


0.88 ±0.02 


991216 


0.25 


~ 10 


1.25 ±0.16 


1.65 ± 0.12 


0.57 ±0.08 


1.57 ±0.14 






0.84 ±0.02 


000301 


0.3 


10 


1.00 ±0.13 


2.83 ± 0.12 


0.99 ±0.05 










000926 


0.4 


10 


1.57 ±0.11 


2.21 ± 0.06 


1.31 ±0.06 




2.06 ± 0.22 


0.96 ± 0.13< c ) 
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1.34 ±0.05 
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0.85 ±0.04^) 


0.63 ±0.04 
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0.2 


8 


0.89 ±0.03 


2.37 ±0.06 


0.93 ±0.05 




3.4 ± 1.5 


1.04 ±0.20( 9 ) 


1.02 ±0.05 



(1) ratio of epochs of first optical measurement and light-curve break ; (2) ratio of epochs of break and last accurate measurement 
of optical transient brightness ; (3) pre-brcak optical light-curve index (F„ oc t~ ai ) ; (4) post-break optical light-curve index ; 
(5) observed slope of optical SED (F v oc v~P) ; (6) pre-break X-ray light-curve index ; (7) post-break X-ray light-curve index (f 
relying on a single, 10 day measurement) ; (8) slope of X-ray SED (Rcfs: Nicastro et al. (1999), Kuulkcrs ct al. (2000), ' c ' 
Harrison et al. (2001) and Piro ct al. (2001), < d > in't Zand ct al. (2001), < e ) Reeves et al. (2003), W Butler ct al. (2003), < 9 > Klose 
et al. (2004) ; (9) optical to X-ray SED slope (eq.[28]) 
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Figure 2. Histogram of the break magnitude Aa, the increase of 
the afterglow power-law decay index, for the optical light-curves 
of the 10 afterglows considered in this work. At the top are in- 
dicated the values expected theoretically for structured outflows, 
jets, and refreshed shocks (energy injection model), interacting 
with homogeneous media (s = 0) and free winds (s = 2), for a 
cooling frequency (u c ) below or above the optical domain (u ). 
Vertical bars indicate upper limits. For structured outflows they 
correspond to a direction toward the observer within the open- 
ing of the uniform outflow core; otherwise Aa can be larger. For 
energy injection and s = 2, the upper limits are for an electron 
energy distribution (cq. [2]) with p = 3.5, while the s = medium 
covers the entire Aa range shown. For jets, the ranges of Aa 
shown correspond to p G (2.0, 3.5). 



the properties of the optical emission (qi, a 2 , and /3 of Ta- 
ble 1) of the GRB afterglows shown in Fig. 1. To take into 
account the uncertainty of the measured light-curve indices 
and SED slopes, we set the following requirements for ac- 
ceptable models. 

The observed SED slope, f3, should be larger or equal 
than the intrinsic one, /3 - Using equations (4) and (6), we 
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Figure 3. Optical slope P for the afterglow 980519, at vari- 
ous epochs (in days). Icr uncertainties of the slope are given 
in parentheses. Fluxes are corrected for Galactic extinction of 
E(B — V) = 0.27. Solid lines are the best power-law fits; the dot- 
ted lines show the power-laws of slopes Icr above and below that 
of the best fit. 



calculate the model SED slope /3 mo dei from the pre-break 
light-curve index «i for models SO, J/O and J, and from 
the post-break index 02 for model EI, and require that 



constraint 1 : (3 > f3 modd - 3 o((3 - f3 mode i) , 



(29) 



where a is the uncertainty of the slope difference. 

For the SO model, the measured break magnitude Aa 
should be smaller than the maximum allowed Aa mo£ j e j given 
in equation (eq.[16]), which is equivalent to requiring that 
q < q (eq.[15]): 

constraint 2a : Aa < Aa mode i + 3 a(Aa - Aa mode i) , (30) 

Given that, for jets, the transition to the post-break de- 
cay takes almost a decade in time for a homogeneous 
medium and two decades for a wind-like medium (Kumar 
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Figure 4. Afterglow 990123, E(B ~V) = 0.02 







0.62d, 


3=0.38(.20) i> 


• 




'o 




0.90d, 


P=0.60(.18) r 










1.60d, [ 


=0.76(.21) h 






E 
























'o 




3.85d, 


3=0.71 (.22) fi 
















• 




- 990510 




I 


R V 


B 



10' 



10 



v(Hz) 



Figure 5. Afterglow 990510, E(B - V) = 0.20 
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Figure 6. Afterglow 991216, E(B - V) = 0.63 
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Figure 7. Afterglow 000301, E(B - V) = 0.05 
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Figure 8. Afterglow 000926, E(B - V) = 0.02 
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Figure 9. Afterglow 010222, E(B - V) = 0.02 
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Figure 10. Afterglow 011211, E(B - V) = 0.04 
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Figure 11. Afterglow 020813, E(B - V) = 0.10 
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Figure 12. Afterglow 030226, E(B - V) = 0.02 



& Panaitescu 2000) and that, for most afterglows, the post- 
break light-curve is monitored for about one decade in time, 
it is possible that the measured decay index «2 underesti- 
mates its asymptotic value. For this reason, we impose the 
following constraint for the J/O and J models: 



constraint 2b : Aa„ 



-5 fT a < Aa < Aa motfci +3 <ja ,(31) 



where a a = a(Aa - Aa mode; ) and Aa mode! is that given by 
equations (eq.[16]) and (eq.[17]), respectively. 
For the EI model, the observed Aa determines the injection 
parameter e (eq.[10]), for which the only constraint is e < 
3 - s. 

A third condition is that the exponent p of the power- 
law electron distribution with energy (eq.[2]) required by the 
observed ai for the SO, J/O, and J models, and by the Q2 
for the EI model should be larger than 2 



constraint 3 : p + 3 a(p) > 2 



(32) 



This condition is set to avoid an electron energy distribu- 
tion with a cut-off or bend at higher energies (otherwise the 
total electron energy diverges). Such an electron distribu- 
tion would provide an explanation for why the radio decay 
of some afterglows is significantly slower than at optical fre- 
quencies (Li & Chevalier 2001, Panaitescu & Kumar 2002), 
but here we choose to ignore this somewhat ad-hoc scenario. 

The models which satisfy the constraints 1, 2, and 3 
above are given in Table 2. Further testing of the viable 
models is done using the X-ray light-curve decay index, S, 
the optical to X-ray slope, f3 ox , and the X-ray SED slope, 
/3 X . If the optical and X-ray decay indices are equal, then 
the cooling frequency u c is not between the optical and X- 
ray domains (with the exception of the J model at t > tb, 
in which case v c is constant in time). Consequently, the X- 
ray flux should lie on the extrapolation of the optical SED 
to X-rays and the slopes of the optical and X-ray continua 
should be the same: 



constraint 4 : S — a =>- j3 — /3 ox = /3 X 



(33) 



If the optical and X-ray decay indices are different, then u c 
is between the two domains, in which case the X-ray flux lies 
below the extrapolation of the optical SED and the X-ray 
spectral slope is larger than that in the optical (eq.[6]): 

constraint 5 : 8 a /3 a < (3ox < /3 X and (3 X = /3 D +l/2 .(34) 

Furthermore, in this case equation (4) leads to 



constraint 6a : <5 



1/4 
-1/4 



s = 
s = 2 



(35) 



for the EI model at t > tb, the SO and J models at t < tb, 
and the J/O model at all times. For the EI model at t < tb 



constraint 6b : S — a = 
while for the SO model at t > tb 
constraint 6c : S — a 



(l + e)/4 
-(l + e)/4 



s = 
s = 2 



(2-q)/(8-q) , 
-(2 - q)/(8 - 2q), 



(36) 



(37) 



Equations (35) - (37) assume that the electron radiative 
cooling is dominated by synchrotron losses. If electrons cool 
mostly through inverse Compton scatterings, a situation 
which we ignore in this work, then S — a is smaller. 



© 2005 RAS, MNRAS 000, 000-000 



Jets, Structured Outflows, and Energy Injection in GRB Afterglows 9 



Table 2. Possible models for the 10 GRB afterglows of table 1 



GRB 


Model 


s 


e or q 


V 


Ho 


Ay 


Note 


Requirement from X-ray 






(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 




980519 


[EI] 





e ~ 


U.O 


^ 3.7 


Z 1.4 





(3 < (3 at 2.1 a 


IC in X-ray 








2 


e ~ 


0.5 


<j O.l 


iC l.U 







inconclusive 






(SO),(J/0),[J] 





g ~ 


0.8 


3.5 


1.25 





(3 < (3 at 1.2 a 


v x < v c 






so,j/o,j 


2 


<j ~ 


1 A 
1.4 


2.8 


0.92 


0.1 




v c < v x 




990123 


EI 


2 


e ~ 


0.5 


2.5 


0.76 







Vx < v c 






[SO],[J/0],{J} 





g ~ 


0.8 


2.6 


0.78 





13 < Po at 2.1 a 


v c < v x 






so,j/o,j 


2 


q ~ 


1.6 


1.9 


0.45 


0.1 




8 l -a l > -1/4 ICW (3 ox 


> (Bo 


99051 


J 


o 






2.2 


0.61 











991216 


[EI] 


2 




0.5 


2.5 


0.77 





(3 < f3 at 2.4 cr 


Vx < Vc 






[SOI [J/Ol u\ 







0.7 


2.7 


0.84 





13 < (3 at 2.1 a 


*'C ■ " X 






SO T /O ( 1\ 


2 




1.4 


2.0 


0.50 





T- 1 An ^ rt, 4- 1 at 1 5 rr 


— a-, ~> —1 /d TPW 


^ Po 


000301 


[J] 









2.3 


0.67(°> 


0.1 


3 Aa > Qi + 3 at 2.4 <r 


no X-ray data 




000926 


[EI] 







0.4 


3.9 


_i Ar -,(h\ 

1.47 w 





/3 < (3 a at 2.2 0- 


IC in X-ray 






EI 


2 


6 — 




o.o 


1.14(b) 


u.uo 




IC in X-ray 






tin t /n / t\ 


u 


<? = 


1.0 


O. 1 


1 flzl(fc) 


n 1 
u. 1 


J. o i_iu <^ ct i -ho at 4 • u 


v x < v c 






Qfi T /O f T^ 

oW, J / W, ^ J ) 


9 


<? = 


1.3 


9 A 


n 79(b) 

U. / Z v ' 


n IK 
u. 10 


J. O LWOi OL\ -7-I tit ± . Z CT 


V c < Vx 






SOf,J/Ot,[J]t 


0,2 


9 = 


1.2 


2.8 


1.38( 6 ) 





J: 3 Aa < ai + 1 at 3.0 cr 


IC in X-ray 




m 0999 




o 


e = 


0.6 


3.4 


1.19 







IC in X-ray 






EI 


2 


e = 


0.9 


2.7 


0.85 


0.1 




IC in X-ray 






(SO),(J/0),{J} 





<? = 


1.3 


2.2 


0.60 


0.2 


SO, J/O: Aa > 3/4 at 1.6 a 


V c < Vx 




011211 


(J) 









2.1 


0.53< c ) 


0.2 


3 Aa > oi + 3 at 1.7 a 


Si-ai> 1/4 at 3.8 a => IC ? 


020813 


EI 





e = 


0.4 


2.8 


0.91 


0.1 




IC in X-ray 






EI 


2 


e = 


0.7 


2.2 


0.57 


0.25 




~ Vx 






(EI)t 


0,2 


e = 


0.5 


2.5 


1.24 





13 < (B a at 1.4 ct 


IC in X-ray 






SO,{J/0} 





9 = 


1.0 


2.1 


0.53 


0.30 


J/O: Aa < 3/4 at 3.0 a 


v c < v x 




030226 


[J] 









2.2 


0.59 


0.1 


3 Aa > «i + 3 at 2.5 ct 


v c < v x 





(1) EI = energy injection, SO = structured outflow, J/O = jet confined by outer outflow J = jet expanding sideways; for unbrackctcd 
models the constraints 1 and 2 are satisfied within 1 <r; for models shown in ( ) and [ ] parentheses, the requirement given in the column 
(7) is satisfied at the la — 2a and 2 a — 3 a levels, respectively; { } parentheses indicate J/O and J models for which the measured 02 — ai 
is 3 a — 5 a below the theoretically expected value ; (2) s = 0: uniform circumburst medium, s = 2: wind-like medium ; (3) parameter 
for energy injection (eq. [10]) or outflow structure (eqs.[13] and [15]), inferred from equations (23), and (25) ; (4) power-law exponent of 
electron distribution with energy (eq.[2]), resulting from equation (4) (for 980519, the <; symbol is followed by p — a p ) ; (5) slope of the 
intrinsic optical spectral energy distribution; for comparison, the following intrinsic afterglow SED slopes and host extinctions were inferred 
from power-law fits to the optical SED and a host SMC-likc reddening curve: (<0 (3 a = 0.57 ± 0.02 (at 3 days) and Ay = 0.14 ± 0.01 
(Jensen et al. 2001), W (3 a = 1.00 ± 0.18 (at 1 day) and Ay = 0.18 ± 0.06 (Fynbo ct al. 2001), < c > (3 = 0.56 ± 0.19 (at 0.6 days) and 
Ay = 0.08 ± 0.08 (Jakobsson ct al. 2003) ; (6) dust extinction in host frame inferred from equations (6) and (9) ; (8) the optical emission 
properties are incompatible with those at X-rays (constraints 4—6) for the models in bold face (ICW = "in contradiction with") . f the 
cooling frequency is below (rcdward of) the optical domain and the electron radiative cooling is assumed to be weaker than synchrotron 
losses; for all other models, the cooling frequency is above (blucward of) the optical. 



If the X-ray afterglow emission is mostly due to in- 
verse Compton scatterings, then equations (33) - (37) are 
not valid. Such a case is easily identifiable through that the 
X-ray emission is above the extrapolation of the optical SED 

Pox < O => Inverse Compton (IC) in X — ray . 

However, a substantial contribution to the X-ray emission 
from inverse Compton scatterings cannot be excluded when 
Pox > Po because the cooling frequency may lie below the 
X-ray domain. 

Whenever p ox > Po, we impose the conditions given in 
equations (33) - (37) at the 3 a level (just as for the optical 
constraints). As shown in Table 2, only a few models are 
eliminated by the X-ray constraints. 



6 CONCLUSIONS 

The first conclusion to be drawn from Table 2 is that an 
energy injection (the EI model) , mediated by some initially 
slower ejecta which catch up with the shock-front (the "re- 
freshed shock" model), may be at work in many afterglows, 
although not in all the afterglows with breaks. For the same 
electron index p, a wind-like medium yields a faster light- 
curve decay than a homogeneous one, which makes the free 
wind medium more often compatible with the steepness 
of the post-break decay. The break magnitude Aa deter- 
mines the exponent of the assumed power-law energy injec- 
tion law (eq.[10]). For a homogeneous medium, we obtain 
0.3 < e < 0.6, while for a free wind, 0.5 < e < 0.9, which 
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imply that the total energy output until the break time (h) 
is larger than the post-burst energy of the GRB ejecta by a 
factor of at least (i(,/ti) e — 1.5 — 4 (2 — 10 for a free wind), 
where t\ is the epoch of the first optical measurement (see 
2nd column of Table 1). If the afterglow energy after the 
end of the energy injection process is close to that previ- 
ously inferred by us from afterglow fits using the jet model 
(Panaitescu & Kumar 2002) , then the GRB output is about 
80 per cent of the ejecta initial energy. Such a high efficiency 
of the GRB mechanism is a serious problem (Kumar 1999) 
for the internal shock model. 

The resulting range for the e parameter and equation 
(11) imply that the distribution of energy with Lorentz fac- 
tor is dE/dT oc r~ 7 with 7 £ (2,3) for a homogeneous 
medium and 7 > 5 for a free wind. The numerical calcula- 
tions of Aloy et al. (2000) for the jet propagation through a 
collapsing massive star yield at break-out dE/dT oc e~°' 7F 
for T £ (5,30). A log-log fit to the log(dM/dr) - T plot 
shown in their Fig. 4 yields dE/dT oc r~ 5 ' 3 for T e (5, 15) 
and dE/dT oc T' 11 for T e (10,30). Therefore the distri- 
bution of ejecta energy with Lorentz factor required by the 
EI model is compatible with that predicted by Aloy et al. 
(2000) only if the circumburst medium has a wind-like strat- 
ification. 

In contrast to the EI model, where the break magni- 
tude Aa is determined by the injection exponent e and, 
hence, can have an arbitrary distribution, a structured out- 
flow (extending this concept to the extreme case of a jet 
with sharp edges) can yield a bimodal distribution of Aa: 
Aa < 3/4 arises when the outflow envelope outshines the 
core (the SO model), or when it is important only dynam- 
ically, by impeding the core's lateral spreading (the J/O 
model), while afterglows with Aa > 3/4 arise from outflows 
with a sufficiently tenuous envelope that the core can expand 
sideways (the J model). The current sample of 10 afterglows 
with light-curve breaks is too small to test the existence of 
a bimodal Aa distribution. If the tentative gap at Aa ~ 1 
shown in Fig. 2 is real, it will favour a homogeneous medium 
around GRBs. 

The second conclusion is that the SO model may be at 
work in more than half of the analyzed afterglows. In these 
cases, a dE/dd oc fl- 1 - * - 3 angular distribution of the ejecta 
kinetic energy per solid angle is required for a homogeneous 
medium and dE/dQ oc _1 - 3 ± 01 f or a f ree w hid. The nu- 
merical simulation of jets produced in type II collapsars, 
shown in fig. 11 of MacFadyen et al. (2001), shows that, af- 
ter having propagated through the helium core, the angular 
distribution of the jet kinetic energy can be approximated 
as dE/dQ oc 6~ A for 6 G (5°, 15°). Such a steep fall-off of the 
ejecta kinetic energy away from the outflow axis corresponds 
to the J/O and J models. 

The third conclusion is that a spreading jet interacting 
with a homogeneous medium is the only viable model for 
the afterglows 990510, 000301, 011211, and 030226, as their 
break magnitudes Aa are too large to be accommodated by 
a structured outflow (assuming the observer to be within the 
opening of the core) or by a jet interacting with a free wind. 
In fact, for the afterglows 000301 and 030226, the measured 
Aa is so large, that even the J model is only marginally 
consistent with the observations. Larger break magnitudes 
would result from combining the EI and J models, though 
this would require a coincidence: the simultaneity of the en- 



ergy injection cessation with the jet edge becoming visible 
to the observer. Stronger light-curve breaks would also re- 
sult from structured outflows if the direction toward the ob- 
server were outside the core opening (Panaitescu & Kumar 
2003), a case not included in this work. Such locations could 
also yield a cusp in the light-curve when the more energetic 
core becomes visible, i.e. at the time of the light-curve break 
(Granot & Kumar 2003). 

A more accurate comparison between the features of 
the above four models for light-curve breaks and the after- 
glow temporal/spectral properties requires numerical cal- 
culations of the outflow dynamics and emission. Such cal- 
culations can follow the transition between the asymptotic 
power-law light-curve decays and may rule out some of the 
models listed in Table 2, if this transition is too slow com- 
pared with the observations. Furthermore, numerical calcu- 
lations provide the framework for inclusion of inverse Comp- 
ton scatterings, which may affect the time-dependence of the 
cooling frequency (and, hence, the light-curve decay index 
at frequencies above it) and may contribute to the X-ray af- 
terglow emission. These two factors have been left out in the 
derivation of the results presented in §2, §3, and §5. Finally, 
numerical modelling of the multiwavelength afterglow emis- 
sion can determine other afterglow parameters (ejecta en- 
ergy, circumburst density, shock microphysical parameters) 
in addition to testing the various possible mechanism for 
light-curve breaks. In this work we have described a simpler 
way to test those mechanisms and have shown how some out- 
flow properties (e.g. energy distribution with Lorentz factor 
or with polar angle) are constrained by afterglow observa- 
tions. 
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